OBSTRUCTIONS TO DETERMINANTAL REPRESENTABILITY 



PETTER BRANDEN 

Abstract. There has recently been ample interest in the question of which 
sets can be represented by linear matrix inequalities (LMIs). A necessary con- 
dition is that the set is rigidly convex, and it has been conjectured that rigid 
convexity is also sufficient. To this end Helton and Vinnikov conjectured that 
any real zero polynomial admits a determinantal representation with symmet- 
ric matrices. We disprove this conjecture. By relating the question of finding 
LMI representations to the problem of determining whether a polymatroid is 
representable over the complex numbers, we find a real zero polynomial such 
that no power of it admits a determinantal representation. The proof uses re- 
cent results of Wagner and Wei on matroids with the half-plane property, and 
the polymatroids associated to hyperbolic polynomials introduced by Gurvits. 



1. Representing sets with linear matrix inequalities 

Motivated by powerful techniques commonly used in control theory, there has 
recently been considerable interest in the following question. 

Question 1. Which subsets of M" can be represented by linear matrix inequalities 
(LMIs)? That is, which sets y are of the form 

y — {(xi, . . . , Xn) e K" : Ao + Aixi + • • • + AnXn is positive semidefinite}, (1) 

where Aq, . . . , A„ are real symmetric m x m matrices? 

In two variables such sets were characterized by Helton and Vinnikov [71 by 
so called rigidly convex sets, thereby answering a question posed by Parrilo and 
Sturmfels [T3]. We will always assume that is in the interior of y and then the 
existence of a LMI representation of y is equivalent to the existence of a monic 
LMI representation, i.e., a representation in which Aq is the identity matrix. 

A polynomial p € M.[xi, . . . , x„] is a real zero polynomial (RZ polynomial) if for 
each X £ M" and p £ C 

p{px) = implies fi is real. (2) 

A set y C M" is rigidly convex (at the origin) if there is a RZ polynomial p for 
which y is equal to the closure of the connected component of 

{a: e M" : p{x) > 0} 

containing the origin. 
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In what follows / will always denote the identity matrix of appropriate size. It 
is not hard to see that if Ai , . . . , A„ are symmetric or hermitian matrices of the 
same size then the polynomial det(/ + xiAi + • • • + x„yl„) is a RZ polynomial. In 
two variables Helton and Vinnikov provided a converse to this fact. 

Theorem 1.1 (Helton-Vinnikov, ^7j). Let p[x,y) be a RZ polynomial of degree d, 
and suppose that p(0, 0) = 1. Then there are symmetric matrices A and B of size 
d X d such that 

p{x, y) = det(/ + xA + yB). 

Theorem 1 1 . 1 1 also settles a conjecture of Peter Lax which asserts that any hyper- 
bolic degree d polynomial in three variables can be represented by a determinant, 
see [T^. By a simple count of parameters one sees that the exact analog of The- 
orem |1.1| in three or more variables fails. However, the count of parameters does 
not preclude a determinantal representation with matrices of a size larger than the 
degree. To this end, Helton and Vinnikov [7, p. 668] made the following conjecture. 

Conjecture 1.2 (Helton-Vinnikov). Let p{xi, . . . , Xn) be a RZ polynomial, and 
suppose that p(0) = 1. Then there are symmetric matrices Ai, . . . , An such that 

p{xi, . . . ,Xn) ^ dct(/ + XiAi H h XnAn). 

In Section |2] we will find a family of counterexamples to Conjecture |1.2[ The 
following relaxation of Conjecture |1.2| has also been suggested. 

Conjecture 1.3. Letp{xi, . . . ,Xn) be a RZ polynomial, and suppose that p(0) = 1. 
Then there are symmetric matrices Ai, . . . , An, and a positive integer N such that 

p{xi, . . . , Xn)'^ = det(/ -I- xiAi H h a;„A„). 

In Section [3] we find a counterexample to Conjecture |1 . 3| by relating the problem 
of finding determinantal representations to the problem of determining whether 
a given polymatroid is representable (comes from a subspace arrangement). The 
counterexample arises from the fact that the Vamos cube is a matroid that has the 
so called half-plane property but is not representable over any field, see [16 . 

The conjecture that any rigidly convex set can be represented by LMIs still 
remains open. 

2. Counterexamples by a count of parameters 

A homogeneous polynomial h{xi, . . . , Xn) G . . . , x„] is hyperbolic with re- 

spect to e e M" if h{e) ^ and if for each x E M" and /i e C 

h(x + /ie) = implies /i is real, (3) 

see [51 [15]. Clearly, if h is hyperbolic with respect to e, then h{x + e) is a RZ 
polynomial. The hyperbolicity cone of /i at e is the set of all x G M" for which the 
univariate polynomial 1 1—^ p{x + te) has only negative zeros. 

We will use the Cauchy-Binet theorem. Let [n] = {1, . . . , n} and let ('^'') denote 
the set of all fc-element subsets of [n] . If A is an n x m matrix and S C [n], T C [m] 
are two sets of the same size we denote by ^4(5*, T) the minor of A with rows indexed 
by S and columns indexed by T. 
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Theorem 2.1 (Cauchy-Binet). Let A be an m x n matrix and B an nxm matrix. 
Then 

det{AB)= A{[m],S)B{S,[m]). 

Theorem 2.2. Let h{x) € ]R[xi, . . . , Xn] he a hyperbolic polynomial with respect to 
e, and let p{x) he the RZ-polynomial defined by p{x) = h{x + e). Lf p admits a 
representation 

p{x) = det(/ + xiAi H + x„An) 

where Aj is symmetric (hermitian) and of size N x N for all j, then p admits a 
representation 

p{x) = det(7 + xiBi H h a;„B„) 

where Bj is symmetric (hermitian) and of size d x d for all j, and d is the degree 
of h and p. 

Proof. By considering a linear change of variables we may, and will, assume that 
h{x) has hyperbolicity cone containing M^^, where M+ is the set of positive reals, 
and that e = (1, . . . , 1)^ =: 1. 

We claim that Aj and I — J2^j^i positive semidefinite (PSD) for all j. The 

univariate polynomial 

det{I + tAj) =p(0,...,f, ...,0) =h{l + tSj), 

where dj is the jth standard bases vector, has only are real and non-positive zeros 
(since 6j is in the closure of the hyperbolicity cone of h). Hence Aj is PSD. Similarly 

det ^I + I-J2A}j = {l+tfp{-l/{l+t),...,-l/{l+t)) 

= (1 + i)^/i(l- 1/(1 + 1/(1 + 

= {l + t)^-'^t'^. 

Hence I — Y^"^i Aj is PSD of rank N — d. Suppose that Ai has rank . Write Ai as 
Ai = X^jLi A-ij, where Aij = VijV*j is PSD of rank 1, v*j is the hermitian adjoint of 

t'ij, and Vij S C^. Similarly we may write / — Aj as a sum X^^^'' Cj, where 

Cj = UjU* is PSD of rank 1. Let now 

(N-d 
Cjy,+Y.A,,x 

where x = {xij)ij and y = (j/i, . . . , j/jv-d) are new variables. Let B be the ma- 
trix with columns u\,. . . ,UN-d,v\i, . . . ,wiri, • • •t'nr„- Rename the columns and 
variables so that B = [wi,...,wm], and the corresponding variables are z = 
{zi, . . . , Zm)- By construction and the Cauchy-Binct theorem 

P{z) = dei{BZB*)= J2 \BiS)\^l[zj, 

sg('5J0 

where Z = diag(zi, . . . , zm), and B{S) = B{[N],S) is the TV x iV minor of B with 
columns indexed by S. 
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We obtain h{x) from P{z) by setting yj ~ 1 and Xij — Xi for all i and j. Since 
all coefficients of P{z) are nonnegative and h{x) is homogeneous of degree d we 
have that for each S with B{S) 7^ there are precisely d indices that correspond to 
X- variables and N — d variables that correspond to ?/- variables. This means that Ur\ 
V = (0), where U — spanjiti}^^'^ and V = spanjuy : \ < i < n and 1 < j < r-i}. 
Hence we may write B as B = PM where P is invertible and M is a block matrix 

[ M2 

d — N columns. Thus 



where Mi has N — d columns and A/2 has M 
P{z) = dct{PP*)dct{MZM*) 



= det(PP* 



E 



E 



iA/i(5i)nM2(52)p n n 



Sie(>j;iJ)5.G(i"-''+^'-i) 
dct(PP*)|det(Mi)|2yi 



jeSi j£S2 



VN-d ^ 



iM2(52)p n 



= det(PP*)| det(Afi)pyi • • • yA,_rfdet 

where mi is the ith column of A'/2- Setting yi 
obtain a representation 



fM-n+d 

ZN-d+imifn* 

\ 1=1 



1 and Xij = Xi for all i and j we 



hix) = det(PP*)|det(A//i)|2det 



where each Ti is PSD of size d x d, and 51]"= 1 positive definite. It follows that 
p{x) has a representation of the desired form. □ 

Nuij [12 proved that the space of all hyperbolic polynomials of degree d that are 
hyperbolic with respect to e G M" has nonempty interior. Hence so does the space of 
RZ polynomials considered in Theorem 2.2 Since any such polynomial that admits 
a determinantal representation also admits a determinantal representation with 
matrices of size d, a count of parameters provides counterexamples to Conjecture 



3. Representability of polymatroids 

We will see here that Question [T] is closely related to the old problem of deter- 
mining if a polymatroid is represent able over C. 

An (integral) polymatroid on a finite set i? is a function r : 2^ ^ N such that 

(1) r(0) = 0; 

(2) If 5 C T C E, then r{S) < r{T)- 

(3) r is submodular, that is, 

r{S UT) + r{S DT) < r{S) + r{T), 
for all subsets S and T of E. 
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A natural class of polymatroids arises from subspace arrangements. Let i? be a 
finite set and V = {Vj)j^E ^ collection of subspaces of a finite dimensional vector 
space V over a field K. Then the function ry : 2^ — > N defined by 



where J^ies ^« smallest subspace containing Ui^sVi, is a polymatroid. This 

follows from the dimension formula for subspaces 



We say that a polymatroid, r, is representable over the field K if there is a subspace 
arrangement V of subspaces of a vector space over K such that r = ry- There are 
several inequalities known to hold for representable polymatroids, see [UlSlin]. The 
simplest of these are known as the Ingleton inequalities. 

Lemma 3.1 (Ingleton inequalities |8]). Suppose that V = (Vi, . . . , Vn) is a subspace 
arrangement. Then 



rv(5i U S2) + rv{Sx U 5*3 U Si) + rviS^) + rv{Si) + rv(S'2 U 5*3 U ^4) < 
rv(5i U S-g) + rv(S'i U S'4) + rv{S2 U ^3) + rv(S'2 U S^) + rv(-5'3 U ^4) 
for all Si, S2, S^, Si C [n] . 



To see that subspace arrangements over C or E are closely related to determi- 
nantal representability we proceed to express the rank function in terms of deter- 
minants. Suppose that Ai, . . . ,A„ are positive semidefinite matrices of the same 
size m, and let V = (Vi, . . . , Vn) be the subspace arrangement in C™ defined by 
letting Vi be the image of Ai for all i. Then 



for all S C [n]. To see this it is enough (by spectral decomposition) to consider the 
case when all matrices are of rank one and that S = [n]. Write Ai as Ai = ViV* 
where Vi £ C™, and let D be the {m + n) x (m + n) diagonal matrix with the first 
m entries equal to one and the remaining entries equal to t. Let further B be the 
m X (m + n) matrix with columns Si, . . . , 6rn,vi, . . . ,Vn, where 6i is the ith standard 
bases vector of C™ . Then by the Cauchy-Binet theorem 



Hence the degree of the above polynomial is the size of a maximal linearly inde- 
pendent subset of {vi, . . . , Vn}, that is, the dimension of Vi + • • • + Vi. 

Next we will see how polymatroids arise from hyperbolic polynomials. This 
connection was observed by Gurvits [6 . If h{xi, . . . ,a;„) is a hyperbolic polynomial 
with respect to e, we define a rank function rankh : M" — > N by 



The rank does not depend on the choice e, but only on the hyperbolicity cone of 
h, that is, deg(/i(e -I- xt)) = deg{h(e' + xt)) for all e' in the hyperbolicity cone 




dim([/ + W)+ dim(J7 nW)^ dim(C/) + dmi{W). 




dctil + tY,A]^det{BDB*)= ^ |B(S')ptl^n{ 
\ i / se(i"+"i) 




m+l,...,m+n}| 



rank;i(a;) = dcg(/i(e + xt)). 



containing e, see [B] and Section [3] 
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The next proposition follows from the work of Gurvits [6_ . He uses Theorem [TT] 
In Section [4] we give a proof that does not rely on the Lax conjecture. 

Proposition 3.2. Let h G ]R[a::i, . . . , Xm] be a hyperbolic polynomial with respect to 
e e M™, and let £ — (ei, . . . , e„) be a tuple of n vectors lying in the closure of the 
hyperbolicity cone of h containing e. Then the function rg : 2["1 — N defined by 



\ies ) 



is a polymatroid. 

A matroid, Ai, may be defined as a polymatroid for which the rank function 
satifies rj^{{i}) < 1 for all i ^ E. Let be a matroid on E. The set of bases of 
M is 

BiM) ^{SCE:\S\^ VMiS) = r^(i?)}. 
It follows from the equivalent definitions of matroids, see |T3], that 

rM{S)^ma.x{\SnB\:BeB{M)}, (4) 

for all S Q E. The bases generating polynomial of Ai is the polynomial in the 
variables {xi)i^E defined by 



Xj. 



sei3(M) jes 

For i Q E, let di G M.^ be defined by Sj{i) — S{i,j), where 5{i,j) is the Kronecker 
delta. By Q 

VMiS) = dcg (^Hm (^1 + tY,SiJ^ , (5) 

for all S CE. 

Let Vg be the Vamos cube, see [T3]. The set of bases of Vg are all subsets 
of size four in Fig. [T] that do not lie in an affine plane. The Vamos cube is not 
representable over any field. However, its bases generating polynomial is hyperbolic 
with hyperbolicity cone containing E^. This follows from the fact that Vg is a so 
called half-plane property matroid (see [3J) which was proved by Wagner and Wei 

m- 

We are now in a position to establish the counterexample to Conjecture |1.3[ 



Theorem 3.3. Let p{x) = hv^iixi + 1, . . . ,xs + 1). Then 

(1) p{x) is a RZ-polynomial; 

(2) There is no positive integer N such that p{x)^ has a determinantal repre- 
sentation. 

Proof. Wagner and Wei [16] proved that hvs{x) is a stable polynomial, that is, 
hY^{x) is non-zero whenever lm{xi) > for all i. Hence, if a; G and y G K^, 
then the polynomial hy^{x + ty) has only real zeros. Thus hvg{x) is hyperbolic 
with hyperbolicity cone containing M^. As previously noted it follows that p{x) ~ 
hvgix + 1) is a RZ polynomial. 

Suppose that there is an integer > for which 

p{x)'^ = det(/ + xiAi + ■■■+ xsAs), 
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Figure 1. The Vamos cube. 



where Ai is hermitian for all i. As in the proof of Theorem |2.2| it follows that 

hv, {x)^ = det(xiBi + • • • + xsBs), 

where Bi is positive semidefinite of size (8A^) x (8iV) for all i. Of course hvg{x)^ 
is also hyperbolic with the same hyperbolicity cone as hyg{x). By ([5|, the rank 
function of hvg{x)^ with respect to £ = . . . , 5%\ satisfies 

r£{S)^deg (^^^ [^ + tY.^5)j^ =Nrv,{S), 

for all 5 C [8]. Hence there is a subspace arrangement V = (Vi, . . . , Vg) for which 

However, it is known that (and thus also Nry^) fails to satisfy Ingleton's in- 
equalities. This is seen by choosing 

5i={5,6}, 52 = {7,8}, 53 = {1,4}, 54 = {2,3}, 

in the Ingleton inequalities. 

□ 



4. Properties of the rank function of a hyperbolic polynomial 

For completeness we give proofs that do not use the Lax conjecture of the prop- 
erties we use about the rank function associated to a hyperbolic polynomial. We 
show that these properties are simple consequences of known concavity properties 
of stable polynomials and discrete convex functions. 

A step from a G Z" to /3 <E Z" is a vector s g Z" of unit length such that 

\a + s- I3\ < \a- l3\, 

where |a| = X]r=i If s is a step from a to /3 we write a A /3. A set C Z" is 
called a jump system if it respects the following axiom. 

(J): If a, f3 £ J' , a P and a + s ^ J^, then there is a step t such that a + s ~> j3 
and a + s + t E J ■ 
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The support, supp(p), of a polynomial p{x) — X^asN" ^ ' ' ' ^Z" the set 
{a G N" : a(a) 7^ 0}. A polynomial p G C[xi, . . . , a;„] is stable ii p(x) 7^ whenever 
lni{xj) > for all j. Let < be the usual product order on Z", i.e., a < (3 if aj < j3j 
for all j. 

Theorem 4.1 ([2 ). The support of a stable polynomial is a jump system. 

Moreover, if all the Taylor coefficients of the stable polynomial p are nonnegative, 
and a, f3 €z supp(p) with a < j3, then 7 G supp(p) for all a < ^ < /3 . 

We need the following simple property of jump systems. 

Lemma 4.2. If J (1 Z" is a finite jump system and a^fi ^ J are maximal (or 
minimal) with respect to <, then \a\ — |/?|. 

Proof. The proof is by contradiction. Let M be the set of maximal elements j3 of 
J , with |/3| — d maximal. Suppose further that /? G M is of minimal L^-distance 
to the set of all maximal (w.r.t. <) elements a with |a| < d. Let a be a maximal 
element that realizes the above distance to (3. 

Clearly aj > Pj for some j. Thus 5j is a step from /3 to a and P + 6j ^ J . By 
(J), /3' = /3 + i5j + s G for some step s from (3 + 5j to a. Since /3 is maximal, the 
nonzero coordinate in s is negative. Now, = so /3' is maximal (w.r.t. <). 
However, — a| < |/3 — aj which is the desired contradiction. □ 

Lemma 4.3. Suppose that h is hyperbolic with respect to e & M" and that ei, . . . , 
lie in the hyperbolicity cone of h, and eo G M". Then the polynomial 



is stable or identically zero. 

Moreover if additionally h[e) > and eo is in the closure of the hyperbolicity 
cone of e, then all Taylor coefficients of p are nonnegative. 

Proof. By Hurwitz' theorem we may assume that ei , . . . , are in the hyperbolicity 
cone containing e. Assume that a G M™ and /3 G M™. Then 



are nonnegative. Clearly q is hyperbolic (or identically zero) with hyperbolicity 
cone containing M'J_, or equivalently, q is homogeneous and stable. It is not hard 
to prove that such polynomials have nonnegative Taylor coefficients, either using 
Renegar derivatives (TS], or as in [TJ [3]. Hence, the Taylor coefficients of p are 



Lemma 4.4. Suppose that h is hyperbolic with respect to e £ M" and that e' lies 
in the hyperbolicity cone containing e. Then 




p(a + il3) l^eo + ^ a^Cj + i ft e^ J j 7^ 0, 

since the hyperbolicity cone is convex, see |5j[T5]. Thus p is stable. 

To prove the last statement we show that all the Taylor coefficients of 

q{xo, Xm) = h{xQeo H h x^e™) 




nonnegative. 



□ 



deg(/i(e + xt)) = deg{h{e' + xt)) 



for all X G M". 
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Proof. The polynomial p{s,t) = h{x 
degree of h be d. By Theorem 4.1 J' - 

deg{h{e + sx)) = deg(s'*p(s"\ 0)) = d-mm{i : (i,0) e J} = d-min{| 
which does not depend on e. 



se + te') is stable by Lemma 4.3 Let the 
supp(p) is a jump system and by Lemma 4.2 

a e J}, 
□ 



a 



Corollary 4.5. Let h be a hyperbolic polynomial with respect to e £ M™, and let 
£ — (ei,...,en) be a tuple of n vectors lying in the closure of the hyperbolicity 
cone of h containing e. Let further J be the support of the stable and homogeneous 
polynomial h{x\e\ + • • • + a:„e„). Then the rank function associated to £ satisfies 



rsiS) = max < ^ 
lies 



Ui : a e J 



for all S Q [n]. 

Proof. Let p(a;i, . . . , Xn) = h{e + XiCi + • • • + Xn^r, 
has nonnegative Taylor coefficients. Hence 



By Lemma 4.3 p is stable and 




OLi : ^ aibi € supp(p) > 
ies J 



Note that the set of maximal elements (w.r.t. <) of supp(p) is equ al to J . Thus 
the inequality r^iS) < maxj^^g^ai : a G ^j} follows from Lemma 
that a £ J and C [n]. Since G supp(p) and < X^ies*^*^* < a we have by 
that 'Ylii^s'^i^i € J . Hence r^iS) > maxj^^g^a 



Theorem 



4.1 



Suppose 
□ 



We may now prove Proposition |3.2 



Proof of Proposition \3.2\ Keep the notation in the proof of Corollary |4.5[ Then 
i7 is a jump system for which all vectors have constant sum. Such jump systems 
are known to coincide with the set of integer points of integral base polyhedra, 
see [11]. Clearly r£ satisfies (1) and (2) of the definition of a polymatroid. The 
submodularity of 



S I— > max < Ui : a £ J' 
lies J 



holds for every constant sum jump system, see [llj. 



□ 
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